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Validation  and  Application  of  Altimetry-Derived 
Upper  Ocean  Thermal  Structure  in  the  Western  North 
Pacific  Ocean  for  Typhoon-Intensity  Forecast 

lam-Fei  Pun,  I-I  Lin,  Chau-Ron  Wu,  Dong-Shan  Ko,  and  W,  Timothy  Liu 


Abstract — This  paper  uses  more  than  501 N)  colocated  and  near* 
coincident  in-situ  profiles  from  the  National  Oceanic  and 
Atmospheric  AdminLstration/Giohal  Temperature  and  Salinity 
Profile  Program  database  spanning  over  the  period  from  2002 
to  2005  to  systematically  validate  the  sateHite-altimetry-derived 
upper  ocean  thermal  structure  in  the  western  North  Pacific 
ocean  as  such  ocean  thermal  structure  information  is  critical  in 
typhoon -intensity  change.  It  is  found  that  this  satellite-derived 
information  is  applicable  in  the  central  and  the  southwestern 
North  Pacific  (covering  122-170°  E,  9-25°  N)  hut  not  in  the  north¬ 
ern  part  (130-170°  E,  25-40°  NK  However,  since  >  of  the 
typhoons  are  found  to  intensify  in  the  central  and  southern 
part,  this  regional  dependence  should  not  pose  a  serious  con¬ 
straint  in  studying  typhoon  intensification.  Further  comparison 
with  the  U5.  Naval  Research  Laboratory's  North  Pacific  Ocean 
Nowcast/Forecast  System  (NPACNFS)  hydrodynamic  ocean  model 
shows  similar  regional  applicability,  hut  NPACNFS  Ls  found  to 
have  a  general  underestimation  in  the  upper  ocean  thermal  struc¬ 
ture  and  causes  a  large  underestimation  of  the  tropical  cyclone 
heat  potential  (TCHP)  by  up  to  60  kj/cm2,  After  validation, 
the  derived  upper  ocean  thermal  profiles  are  used  to  study  the 
intensity  change  of  supertyphoon  Dianmu  (2004).  It  Ls  found  that 
two  upper  ocean  parameters,  Le,,  a  typhoon's  self-induced  cooling 
and  the  during-typhoon  TCHP,  are  the  most  sensitive  parameters 
(with  R?  ~  0.7)  to  the  6-h  intensity  change  of  Dianmu  during  the 
study  period  covering  Oianmo's  rapid  inteasification  to  category  5 
and  its  subsequent  decay  to  category  4.  This  paper  suggests  the 
usefulness  of  satellite* based  upper  ocean  thermal  information  in 
future  research  and  operation  that  is  related  to  typhoon-intensity 
change  in  the  western  North  Pacific. 

hide x  Terms — Altimetry,  typhoon,  upper  ocean  thermal 
structure. 


I.  Introduction 

THE  WESTERN  North  Pacific  is  among  the  world  oceans 
where  most  tropical  cyclones  occur.  Here,  an  average 
of  30  cyclones  (including  tropical  storms,  depressions,  and 
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Fig.  1.  (a)  Flowchart  showing  ibe  steps  thai  arc  involved  in  the  TLM  WNR 
(b)  Schematic  diagram  illustrating  the  TLM_WNP  method  in  retrieving  upper 
ocean  thermal  profile, 

typhoons  from  category  1  to  51)  are  found  each  year.  In  com¬ 
parison,  ^  20  cyclones  are  found  in  the  Eastern  Pacific  and 
the  Atlantic,  while  S  10  cyclones  are  found  in  the  Southern 
Pacific  Ocean  each  year.  In  the  western  North  Pacific,  cyclones 
profoundly  impact  the  half-billion  people  living  near  the  Asian 
coasts.  This  can  be  clearly  demonstrated  in  the  case  of  re¬ 
cent  Supertyphoon  Saomai  {August  2006)  in  which  eastern 
China  was  devastated  with  a  death  loll  of  3(H)  and  a  loss  of 
1,5  billion  USD. 

Although  with  such  impact,  the  current  typhoon- iniensity- 
forecast  skill  in  the  western  North  Pacific  still  has  much  room 
for  improvement.  Unlike  the  situation  in  the  Atlantic  where 
coupled  hurricane-oeean  models  are  used  in  the  operational 
forecast  [  1  ],  here,  most  of  the  operational  models  remain  un¬ 
coupled  and  use  only  the  prescribed  sea-surface  temperature 
(SST)  to  represent  the  ocean's  role  in  the  forecast.  As  pointed 
out  by  Bender  and  Ginis  ( 1  ],  Emanuel  [3),  Emanuel  et  aL  [4], 
and  Shay  et  aL  [32],  this  is  not  realistic  because  cyclones  inter* 
act  not  only  with  the  surface  but  also  with  the  entire  upper  (Kean 
column  (typically  down  to  200  m).  In  the  uncoupled/prescribed- 
SST  forecast,  the  ocean's  dynamic  response  due  to  the  cy¬ 
clone's  wind  forcing  is  not  included;  thus,  the  critical  negative 
feedback  due  to  the  cyclone's  self-induced  cooling  is  not  con¬ 
sidered  and  can  lead  to  considerable  error  in  intensity  forecast 


lThc  Saffir-Simpson  Tropical  Cyclone  Seale  based  on  a  I ’min  average 
sustained  wind  speed  (tropical  depression:  <  34  kn,  tropical  storm:  34-63  kn. 
category  1 :  64-82  kn.  category  2:  83-95  kn.  category  3:  96-1 1 3  kn.  category  4; 
1 14-135  kn.  and  category  5:  >  135  kn). 
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Fig.  2,  (a)  Monthly  climatological  ratio  between  D26  and  D2f>  during  May-Ociober  in  the  western  North  Pacific  (Dam  source:  NOA A/NODC/WOAO l ,  [3T)>, 

Gray  represents  land  or  areas  where  the  SST  is  lower  than  the  26  °C  and  no  ratio  can  be  calculated,  (b)  Monthly  climatological  mixed-layer  depth  from  the 
U.S.  NRL 
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[1J.  PI  14]*  113].  To  indude  the  ocean’s  dynamic  response 
in  the  forecast,  one  important  prerequisite  is  to  operationally 
obtain  the  upper  ocean  thermal  structure  information  [I],  [4], 
[20],  J31],  [35].  However,  how  to  obtain  such  information  over 
the  vast  western  North  Pacific  ocean  is  a  big  problem  since  it  is 
impossible  to  use  in-situ  measurements  to  meet  this  operational 
need.  The  satellite  remote-sensing-based  method  thus  becomes 
vital.  It  has  been  proposed  by  Goni  et  ai  [9]  and  Shay  et  al. 
[32]  that  by  using  satellite  altimetry  and  SST  data  as  input  to  a 
simple  two- layer  reduced-gravity  ocean  model,  one  can  obtain 


a  first -order  approximation  of  the  upper  ocean  thermal  infor¬ 
mation.  This  method  has  the  advantage  of  fast  calculation  and 
simple  implementation;  thus,  it  is  suitable  for  the  operational 
needs  and  is  currently  adopted  by  the  National  Oceanic  and 
Atmospheric  Administration  (NOA A)  to  operationally  provide 
upper  ocean  thermal  information  over  global  (Keans  [10]. 

Clearly,  it  is  necessary  to  validate  this  derived  information  so 
as  to  assess  its  applicability  regionally.  Since  there  is  little  such 
information  for  the  western  North  Pacific  region,  this  paper 
aims  lo  conduct  systematic  validation  using  a  large  set  (>  5000) 
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of  coincident/coloeated  sea  truth  data  from  ihe  NOA  A/Global 
Temperature  and  Salinity  Profile  Program  (GTSPP,  [  17])  during 
2002-2005  to  contribute  to  this  need.  Because  this  is  a  rela¬ 
tively  simple  model,  additional  comparison  with  a  26- sigma 
level  operational  ocean  model  (the  North  Pacific  Ocean 
Nowcast/Forecast  System  (NPACNFS)  from  the  U.S.  Naval 
Research  Laboratory  (NRL)  [181)  is  made.  After  validation,  we 
use  the  derived  upper  ocean  thermal  structure  information  to 
study  the  intensity  change  of  the  most  rapidly  intensified  (from 
category  l  to  5  in  18  h)  western  North  Pacific  typhoon  since 
1960,  i,e„  supenyphoon  Dianmu  (2004). 


11,  Upper  Ocean  Thermal  Structure:  That  Is 
Derived  From  the  Two-Layer  Method  and 
Tropical  Cyclone  Heat  Potential  (TCHP) 


Fig.  1(a)  summarizes  the  steps  that  are  involved  to  estimate 
the  upper  ocean  thermal  structure  information  using  satellite 
altimetry  and  SST  data  in  a  two- layer  reduced -gravity  ocean 
model,  as  proposed  by  Goni  et  al  [9]  and  Shay  et  al.  [32),  In 
the  two- layer  scheme,  the  ocean  is  separated  into  two  layers 
of  fluid.  The  20  °C  isotherm  (which  usually  lies  at  the  center 
of  the  main  thermocline  in  most  tropical  regions  [10],  [32]) 
serves  as  a  border  that  separates  the  upper  and  lower  layers. 
The  distance  from  the  sea  surface  to  the  20  °C  isotherm  (D20)  is 
thus  defined  as  the  upper  ocean  thickness  hy.  As  the  variations 
in  the  depth  of  the  main  thermocline  (here  defined  as  D20  or 
h) )  are  associated  with  the  variations  in  the  sea-surface  height. 
Thus,  h\  (or  D20)  can  he  estimated  from  the  observed  altimetry 
sea-surtace  height  anomaly  (SSHA  or  if)  data  as  follows  [9]: 


=  hi(x,y)  + 


P2 (x.y) 


P2[x*v)-Pi[x*y) 


(i) 


where  h\  is  the  mean  climatological  upper  layer  thickness 
(Le„  the  climatological  D20),  and  pi  and  (n  are  the  density 
of  the  upper  (sea  surface  to  20  °C  isotherm)  and  lower  (20  DC 
isotherm  to  ocean  bottom)  layers,  respectively.  In  this  pa¬ 
per,  the  SSHA  data  are  from  the  let  Propulsion  Laboratory's 
delayed  mode  data  from  the  TOPEX/Poseidon  and  Jason- 1 
satellite  altimeters  [7],  With  this  product,  SSHA  is  based  on 
the  difference  between  the  observed  sea- surface  height  and 
the  six- year  (1993-1998)  mean  of  the  sea-surface  height  data 
from  multiple-altimetry  missions  including  TOPEX/Poseidon, 
ERS-1,  ERS-2,  and  GEOSAT  |24).  h\  is  from  the  high- 
resolution  (25  x  25  km)  climatological  temperature  data, 
Le.,  the  WOA01  data  set  (hased  on  the  in-si tu  data  dur¬ 
ing  1982-2001)  of  the  NOAA/National  Ocean  Data  Center 
(NODC)  [33],  pi  and  p2  are  also  calculated  from  the  WOAOI 
temperature  and  salinity  database  [22t,  [34).  Clearly,  all  pa¬ 
rameters  are  functions  of  location  (x,  y)  and  time  ((),  except 
the  climatological  parameters  (i.e.,  h\  and  p),  which  are  only 
functions  of  location  (x,y),  After  estimation  of  the  initial 
D20,  the  depth  of  the  26  °C  isotherm  (D26)  is  estimated 
based  on  the  climatological  ratio  between  D20  and  D26  for 
each  1/4U  by  1/4°  grid  point  in  the  western  North  Pacific 
(domain:  120-170°  E,  9-40°  N)  as 


(2) 


Fig.  3,  Luxations  of  the  four  validation  zones  in  the  western  North  Pacific 
including  the  SNFG.  die  SERZ,  ihe  NERZ.  and  the  1ERZ.  Circles  denote  the 
locations  of  the  5684  GTSPP  in-siiu  profiles  that  were  used  in  the  validation. 


TABLE  t 

Number  or  monthly  In-Shv  GTSPP  Prohlls  That  Arl  Usm  in 
validation  in  Dihtrknt  Zonks  During  2002-2005 
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where  //  represents  D26,  and  0  is  the  climatological  ratio 
that  is  also  calculated  from  the  WOAOI  temperature  database 
[33],  From  [10]  and  [32 J,  it  is  important  to  estimate  the  26  X 
isotherm  because  26  °C  represents  a  threshold  temperature 
for  tropical  cyclone  genesis  and  is  commonly  used  in  tropical 
cyclone  research  [6],  [23].  The  0  maps  for  each  month  during 
the  May-October  western  North  Pacific  typhoon  season  are 
shown  in  Fig.  2(a). 

Next,  SST  (Fig.  I)  can  be  obtained  directly  from  the  cloud- 
penetrating  Tropical  Rainfall  Measuring  Mission/Microwave 
Imager  (TRMM/TMI)  and  the  Advanced  Microwave  Sounding 
Radiometer  for  the  Earth  Observing  System  (AMSR-E)  [36), 
Finally,  the  mixed- layer  depth  is  obtained  from  the  monthly 
averaged  climatological  mixed-layer  depth  data  from  the  U.S. 
NRL  [16]  [Fig.  2(b)],  Thus,  the  upper  ocean  depth-temperature 
profile  can  he  produced  operationally  to  a  first  order  at  1/4° 
(~25  km)  spatial  resolution  in  a  daily  basis  [Fig.  1(b)]  [9],  [10], 
[32].  It  is  referred  to  as  TLM_WNP  (Two- Layer  Method  for 
the  Western  North  Pacific)  hereafter  in  this  paper. 

After  obtaining  the  upper  ocean  depth -temperature  profile, 
TCHP  can  be  calculated.  TCHP  was  originally  proposed  by 
Leipper  and  Volgenau  in  1972  as  the  integrated  heat  content 
from  the  depth  of  26  X  isotherm  to  the  surface.  Recent  works 
by  Goni  and  Trinanes  [10]  and  Shay  et  al  [32]  report  that 
as  TCHP  carries  the  information  of  the  upper  ocean  thermal 
structure,  it  can  provide  useful  information  in  studying  the 
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Hg,  4r  Monthly  scalier  plots  of  the  TLM_  WNP-esti mated  D20  ami  the  in -situ  GTSPP  D20  in  the  (a)  SNFG.  (b)  SERZ,  (c)  IKRZ,  and  (d)  NKRZ  during 
2002-2005,  The  matching  pairs  are  separated  into  three  groups  representing  warm  feature  pairs  ISSHA  >  8  cm-  red  circles),  neutral  pairs  (SSHA  between  -S 
and  S  cm.  green  triangles),  and  cold  feature  pairs  (SSHA  <  —8  cm,  blue  squares). 


cyclone's  intensity  change.  They  proposed  that,  after  obtaining 
the  upper  ocean  thermal  structure  from  the  two-layer  method, 
TCHP  Qh  [101. 119],  [32]  can  be  calculated  as 

Qh  (x,y,t)  =  rpPAT(x,y,z,t)AZ  (3) 

where  cp  is  the  capacity  heat  of  the  water  at  constant  pressure, 
which  is  taken  as  4 1 78  J  kg  1  °C“l;  p  is  the  average  density 
of  the  upper  ocean,  which  is  taken  as  1026  kg  ■  m “3;  and  AT  is 
the  difference  between  T(z)  and  26  °C  over  depth  interval  A Z. 
Clearly,  Qp  is  also  a  function  of  location  (i,  y)  and  time  (f). 

OX.  Validation  of  the 

Two- Layer- Method* Estimated  Upper  Ocean 
Thermal  Structure  by  /*-5m/  Profiles 

In  this  paper,  5684  colocated  and  near-coincident  (ac¬ 
quired  on  the  same  day)  in-sit  u  temperature  profiles  from  the 
NOAA/GTSPP  [17]  covering  the  May  ^October  typhoon  season 
from  2<X>2  to  2005  are  used  for  validation.  Due  to  distinct 
hydrographical  settings  in  the  western  North  Pacific  [14],  [27]. 
[30],  [37],  validation  is  done  according  to  the  four  zones  in  the 
western  North  Pacific.  These  four  zones  are  the  Southern  North 
Pacific  Gyre  (SNPG;  122-170°  E,  9-18°  N);  two  known  eddy- 


rich  zones,  i.e.,  the  South  Eddy-Rich  Zone  (SERZ;  122-170°  E, 
1 8-25°  N)  and  the  North  Eddy-Rich  Zone  (NERZ;  140-170*  E, 
30-40°  N)  [14],  [20],  [27],  [30],  [37];  and  lastly,  the  Intermedi¬ 
ate  zone  between  the  two  Eddy-Rich  Zones  (IERZ;  1 30=170°  E, 
25-30°  N  ).  The  locations  of  these  four  zones,  together  with 
the  locations  of  the  in-situ  profiles,  are  shown  in  Fig.  3.  The 
distribution  of  the  location  of  these  profiles  over  these  zones 
during  the  typhoon  season  is  given  in  Table  I.  The  validation  is 
done  for  the  TLM_WNP  D20,  D26,  and  SST.  For  each  month 
at  each  zone,  the  root  mean  square  (rms)  di (Terence  (which 
represents  the  absolute  difference),  the  Percentage  of  RMS  in 
the  Observed  Mean  (PROM,  which  represents  the  difference  of 
the  estimated  values  relative  to  the  observation),  and  bias  are 
calculated. 

A.  Validation  of  the  Depth  of  the  2(f  C  isotherm  (D20) 

D20  is  the  first  parameter  that  is  estimated  from  the  two- layer 
algorithm  (Section  II)  and  is  used  to  derive  the  26  °C  isotherm 
(D26),  It  is  thus  fundamental  to  assess  the  accuracy  of  D20  first 
since  incorrect  D20  will  lead  to  subsequent  incorrect  estimation 
of  D26.  The  monthly  scatter  plots  comparing  the  derived  and 
in-situ  D20  during  2002-2005  arc  given  in  Fig.  4,  In  addition, 
in  Fig.  4,  the  matching  pairs  arc  separated  into  three  groups 
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Fig.  5.  (a)  Monthly  nms  difference,  (b)  PROM,  and  (c)  bias  of  the 
TLM.WNP-e&timated  D20  in  each  zone.  The  dotted  line  in  (b)  depicts  the 
accuracy  threshold  of  25%  PROM. 

representing  positive  SSHA  pairs  (SSHA  >  8  cm.  correspond¬ 
ing  to  large  values  of  D20),  neutral  pairs  (SSHA  between  -8 
and  8  cm),  and  negative  SSHA  pairs  (SSHA  <  -8  cm,  corre¬ 
sponding  to  small  values  of  D20)  [27],  The  monthly  accuracy 
analysis,  in  rms  difference,  bias,  and  PROM,  are  given  in  Fig,  5, 
As  in  May,  there  are  very  few  matching  data  points  (22  pairs) 
in  the  NERZ  that  can  be  used  for  validation  because  the  water 
temperature  was  still  cold  and  mostly  below  20  °C  and  thus 
are  not  discussed.  From  Fig.  4,  one  can  see  that  throughout 
May-October  in  the  SNPG  [Fig,  4(a)]  and  SERZ  [Fig.  4(b)], 
the  derived  TLM_WNP  D20  is  in  good  agreement  with  the 
in-situ  D20,  However,  moving  northward  to  the  IERZ 
[Fig.  4(c)],  clear  increase  in  spreading  is  found.  In  particular, 
ihc  spreading  is  found  to  be  caused  by  overesiimation  in  the 
positive  SSHA  pairs  (red  circles)  and  underestimation  in  the 
negative  SSHA  pairs  (blue  squares).  Further  north  to  the  NERZ 
1  Fig.  4(d)],  a  huge  error  is  seen.  One  can  see  that  the  derived 
D20  can  reach  300  m  in  the  positive  SSHA  region,  with  an 
overestimation  of  —200  m  from  the  in-situ  D20  of  —100  m 


TABLE  II 

Applicability  op  the  Depth  or  20  °c  Isotherm  (D20)  Using  the 
tlm_wnp  Method  for  Different  Months  and  Zones.  Using  the 

25%  PROM  AS  THE  ACCEPTANCE  THRESHOLD,  THJ  “APPLICABLE” 

Months  and  Zones  are.  shown  by  Ticks,  while  the  Other 
Months  and  Zones  Are  Shown  by  Crosses 
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[also  in  Fig.  6(a)],  while  underestimation  in  t  he  negative  SSHA 
region  is  again  observed. 

From  the  accuracy  analysis  in  Fig,  5,  we  can  see  that,  in  the 
SNPG  and  SERZ,  the  rms  difference  is  usually  about  25  m  or 
less,  bias  <  20  m,  and  PROM  <  25%,  However,  in  the  IERZ, 
the  bias  increases  to  >  25  m,  and  the  PROM  reaches  —30%.  As 
in  the  NERZ,  both  rms  difference  and  bias  can  be  as  much  as 
100  m,  while  the  PROM  can  also  reach  100%.  Taking  a  PROM 
of  25%  as  a  threshold  guide,  the  applicable  regions/months  are 
determined  and  summarized  in  Table  II.  Basically,  the  two- 
layer-  met  hod- derived  D20  is  of  acceptable  accuracy  in  the 
SNPG  and  SERZ  but  not  in  the  IERZ  and,  particularly,  NERZ. 
Possible  reasons  for  these  inapplicabilities  are  discussed  as 
follows:  From  Fig.  6,  one  can  see  that  the  hydrographic  setting 
in  the  NERZ  and  IERZ  is  complex.  There  exists  a  mixed-water 
region  where  the  warm  Kuroshio  current  from  the  south  and  the 
cold  Oyashio  current  from  the  north  mixed,  a  cold  sub- Arctic 
current,  and  a  recirculation  region  in  the  Kuroshio  extension 
[28],  With  this  complex  situation,  the  simple  two-layer  method 
may  not  work.  Another  possibility  is  the  assumption  that  is 
made  in  choosing  D20  as  the  border  separating  the  two  layers 
(Section  II).  D20  was  chosen  and  proposed  originally  as  it 
is  usually  where  the  main  thermocline  in  the  study  region  of 
Shay  et  ai  [32],  i.e.T  the  Gulf  of  Mexico,  lies.  In  1 10],  this 
method  and  the  choice  of  D20  are  subsequently  applied  glob¬ 
ally.  It  is  possible  that,  in  the  IERZ  and  NERZ*  D2Q  may  not  be 
an  appropriate  choice  in  representing  the  thermocline  location 
and  resulted  in  the  observed  error  in  estimation.  As  reported  by 
Qui  and  Chen  [29],  in  these  regions,  the  main  thermocline  is 
below  20  °C  and  can  reach  as  low  as  12  °C,  suggesting  possible 
errors  of  using  D20  in  these  regions.  This  can  be  tested  in  future 
investigations  and  assess  whether  accuracy  can  be  improved  by 
changing  the  choice  of  D20  in  the  IERZ  and  NERZ. 

B  Validation  of  the  Depth  of  the  26  °C  Isotherm  (D26) 
andSST 

After  estimating  D20,  D26  can  then  he  estimated  from  (2) 
(Section  II)  based  on  the  estimated  D20.  According  to  the 
validation  results  of  D2Q  (Section  UFA),  the  two- layer  method 
cannot  produce  reliable  D20  in  the  IERZ  and  NERZ.  Thus, 
in  these  two  zones,  D26  is  not  calculated.  As  in  Fig.  7(a),  it 
can  be  seen  in  the  monthly  scatter  plots  that  the  derived  and 
observed  D26  in  the  SNPG  is  again  in  good  agreement  for  all 
the  months  between  May  and  October.  The  corresponding  rms 
difference  is  typically  about  12  m  with  bias  <  3  m  and  a  PROM 
of  about  12%  (Fig.  8).  In  the  SERZ,  a  similar  agreement  is 
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Fig.  6.  (a)  Observational  temperature  transeel  of  145°  E  in  the  NERZ  (after  f.T7]).  (b)  Hydrographical  setting  of  the  western  North  Pacific  (after  [28]). 


found  from  July  to  October  [Figs.  7(h)  and  8],  In  May  and 
June,  dear  underestimation  in  the  positive  SSHA  pairs  (red 
circles)  is  seen  [Fig.  7(b)].  This  leads  to  a  decrease  in  the 
accuracy  that  the  rms  difference  is  about  22-25  m  and  the 
PROM  is  about  32%-38%.  This  shows  unreliable  estimation 
of  D26  in  these  two  months.  This  inapplicability  of  the  two- 
layer  method  in  the  SERZ  during  the  early  typhoon  season  (e.g.* 
May  and  June)  may  be  caused  by  the  uncertainty  of  0  (i.e„ 
the  climatological  ratio  between  D20  and  D26).  As  shown  in 
Fig.  2(a),  during  these  months,  0  is  of  very  Low  value  (0-0.3) 
and  shoals  sharply  in  the  SERZ  (122-1 70°  E*  1 8-25 :  N)  during 
this  early  period  before  summer.  With  this  very  low  and  sharp 
decrease  in  0 ,  D26  can  be  underestimated  as  observed.  Table  III 
summarizes  the  applicable  months  in  the  SNPG  and  SERZ  of 
the  D26  estimation.  The  validation  results  of  SST  are  given  in 
Fig.  9.  One  can  see  that  the  monthly  rms  difference  between 
the  TLM.WNP  (he.,  the  TRMM/TMI  and  AMSR-E  SST) 
and  the  observed  GTSPP  SST  is  within  the  range  of  0.4  °C- 
I  °C*  The  bias  is  typically  below  0.5  °C,  showing  good  agree¬ 
ments.  A  similar  finding  was  also  reported  by  (36). 

C,  Applicable  Zsmes  and  Periods 

As  summarized  in  Tables  II  and  IIL  the  two-layer  method 
is  applicable  in  the  southern  part  of  the  western  North  Pacific 
ocean  (i.e,,  the  SNPG  and  SERZ)  but  not  in  the  northern 
part  (i.e.,  IERZ  and  NERZ).  Thus*  SNPG  and  SERZ  are  the 
“applicable"  zones  where  this  method  can  be  used  with  ac¬ 
ceptable  accuracy  (typically,  the  monthly  PROM  is  <  25%). 
In  addition,  the  applicable  months  of  the  SNPG  are  from  May 
to  October  (Le.*  the  entire  typhoon  season),  while  those  of 
the  SERZ  are  from  July  to  October.  To  assess  whether  this 
limited  applicability  in  location  and  time  may  pose  a  seri¬ 
ous  constraint  in  studying  typhoon  intensification,  we  analyze 
45  years  (1960-2004)  of  typhoon  track  and  intensity  data  over 
the  western  North  Pacific  during  the  May-October  period.  The 


data  source  of  this  analysis  is  the  best  track  data  of  the  U.S. 
Joint  Typhoon  Warning  Center  (JTWC).  Table  IV  compares 
the  number  of  typhoons  (during  intensification)  passing  the 
applicable  zones  during  the  applicable  months  with  the  total 
(he,,  passing  both  inside  and  outside  of  the  applicable  zones  and 
months).  As  shown  in  Table  IV,  this  limited  applicability  should 
not  pose  a  severe  constraint  in  studying  the  intensification 
of  category  5  (in  Saffir-Simpson  scale)  typhoons  since  100% 
of  category  5  typhoons  are  found  to  intensify  only  in  the 
applicable  zones  [21].  One  can  also  see  that  96.2%  of  category  4 
intensifying  typhoons  and  97%  of  category  3  typhoons  pass 
the  applicable  zones  and  applicable  months.  This  suggests  that* 
although  TLM_WNP  is  only  applicable  in  the  applicable  zones 
during  applicable  months,  this  limited  applicability  is  not  a 
severe  hindrance  when  studying  the  intensification  of  intense 
(i.e.*  categories  3-5)  typhoons  since  most  (96.2%- 100%)  of 
these  intense  typhoons  are  found  to  intensify  in  the  applicable 
zones  during  applicable  months.  The  limited  applicability  of  the 
two- layer  method  for  category  2  and,  particularly*  category  1 
typhoons*  however*  poses  some  constraints.  As  shown  in 
Tabic  IV,  61 ,9%  of  the  category  I  typhoons  pass  the  applicable 
zones  during  applicable  months  during  their  intensification, 
while  the  other  38.1%  of  the  category  l  typhoons  do  not.  As 
for  category  2*  87.9%  of  the  cases  are  applicable.  Fig,  10  sum¬ 
marizes  the  intensification  track  segment  for  all  562  typhoons 
(from  category  1  to  5)  passing  the  western  North  Pacific  and* 
in  average*  84.5%  of  the  typhoons  that  are  intensified  in  the 
applicable  zones  during  applicable  months  (as  also  shown  in 
Table  IV). 

IV,  Comparison  With  thk  NPACNFS 
Full  Ocean  Model 

Since  the  two- layer  reduced-gravity  model  is  a  simple  model* 
here  we  further  compare  iLs  performance  with  the  NPACNFS 
full  ocean  model  from  the  U.S,  NRL  [18].  NPACNFS  is  a  3-D 
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Fig,  7,  As  in  Fig.  4,  the  monthly  scatter  plots  of  the  TLM_WNP-estimated  D26  and  the  in-situ  GTSPP  D26  in  the  (a)  SNPG  and  <b)  SHRZ  during  2002-2005, 


Fig,  H.  (a)  Monthly  rms  difference,  (b)  PROM,  and  (c)  bias  of  the 
TLM_WN P-esli mated  D26  in  SNPG  and  SERZ.  The  dotted  line  in  lb)  depicts 
the  accuracy  threshold  of  259h  PROM. 

hydrodynamic  model  with  26  sigma  levels  in  the  vertical*  The 
model  assimilated  satellite  data.  Multichannel  SSI  (MCSST)* 
and  static  climatology  from  the  Modular  Ocean  Data  Assimi¬ 
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Fig.  9.  Comparison  of  the  value  of  (a)  rms  difference  and  lb)  bias  for  the  SST 
in  the  TLM_WNP  in  the  SNPG,  SERZ,  IERZ,  and  NBKZ. 


lation  System  to  produce  3-D  fields  of  temperature  and  salinity 
[18],  Three  months  (July-September  2003)  of  the  NPACNFS 
upper  ocean  thermal  structure  data  is  obtained  for  compar¬ 
ison.  Again*  the  NPACNFS  data  is  validated  by  the  in- situ 
GTSPP  data. 
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TABLE  IV 

STATISTICS  OF  45  YEARS  (1960-2004)  oi  WESTERN  NORTH  PACIFIC 
Typhoons  (Category  i  Through  5)  Based  on  thi  Best  Track 
Data  From  the:  Joint  Typhoon  Warning  Center.  The  Total 
Numbers  of  Typhoons  (During  Intensification)  are 
Compared  With  thi  Hummers  of  Typhoons  Passing  the 
Applicable  Zones  During  Appijcam  i  Months 
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Fig.  1 0.  Tracks  of  the  562  typhoons  in  the  past  45  years  f  1 960-2004)  passing 
the  western  North  Pacific  and  the  applicable  areas, 

A ,  Depth  of  the  20  °C  Isotherm  (D20) 

The  validation  results  of  the  NPACNFS  D20  are  shown  in 
Fig,  l  l(aHc).  It  can  he  seen  that,  during  the  three  months  of  the 
validation  period,  the  NPACNFS  D20  shows  slight  increase  in 
accuracy  in  Ihe  SNPG  as  compared  with  the  two-layer  method 
(Fig,  5)  that  the  monthly  rms  difference  is  typically  about 
13-20  m  with  a  PROM  <  10%,  A  small  negative  bias  of 
about  10  m  in  the  SNPG  between  July  and  September  is 
noted.  In  the  SERZ,  the  NPACNFS  and  TLM_WNP  have 
similar  accuracies  in  the  D20  estimation,  with  a  PROM  of 
about  I!%-I6%  and  an  17ns  difference  of  about  21  m  in 
August  and  September  [Fig.  11(a)  and  (b)].  A  larger  monthly 
rms  difference  of  32  m  is  noted  in  the  NPACNFS  D20  in 
July  [Fig.  11(a)].  In  the  IERZ,  the  D20  from  the  NPACNFS 
performs  better  than  that  from  the  two-layer  method  (Fig.  5); 
for  example,  in  July  and  September,  the  PROM  is  <  20% 
(rms  difference  <  20  m),  while  in  August,  the  PROM  is  —29% 
[Fig.  11(a)  and  (b)|,  In  the  NERZ,  the  D20  from  the  NPACNFS 
is  better  than  that  from  the  two- layer  method,  with  a  monthly 
PROM  in  the  range  of  43%-56%  and  an  rms  difference  of 
about  30-46  m  [Fig,  11(a)  and  (b)],  as  compared  with  the  two- 
layer  PROM  of  >  90%  and  mis  difference  of  70  m.  However, 
although  NPACNFS  behaves  better  than  the  two- layer  method 


in  the  NERZ,  the  rms  of  —40  m  and  PROM  of  —50%,  as 
validated  by  the  in-situ  points,  are  still  too  high  and  are  still  not 
suitable  for  use.  A  summary  of  the  applicability  for  applicable 
zones  and  months  using  the  NPACNFS  is  given  in  Table  V. 

B .  Depth  of  the  26  °C  Isotherm  ( D26)  and  SST 

The  validation  results  of  the  NPACNFS  D26  are  summarized 
in  Fig.  11{dH0*  1b  comparison  with  the  two-layer  method 
(Fig,  8),  the  NPACNFS  D26  has  little  less  accuracy  in  the  SNPG 
that  the  monthly  PROM  is  >  14%,  as  compared  with  ihe  11% 
of  the  two- layer  method.  In  addition,  a  clear  negative  bias  in 
the  range  of  8%-23%  is  found.  In  the  SERZ,  the  two-layer  D26 
has  higher  accuracy  that  the  PROM  of  the  NPACNFS  D26  can 
be  up  to  46%  [Fig.  1 1(e)],  while  the  two-layer  PROM  is  ahout 
20%  [Fig.  8(b)],  Moving  northward  to  the  IERZ  and  NERZ, 
an  obvious  negative  bias  of  the  NPACNFS  D26  with  PROM 
>  40%  is  found,  again  suggesting  that  in  these  two  zones, 
even  with  the  complex  NPACNFS,  D26  cannot  he  reliably 
estimated.  Again,  taking  the  25%  monthly  PROM  accuracy  as 
a  threshold,  the  applicability  of  the  NPACNFS  is  summarized 
in  Table  VI. 

The  validation  results  of  the  NPACNFS  SST  with  in-situ  data 
are  given  in  Fig.  11(g)  and  (h).  In  contrast  to  the  high  SST 
accuracy  that  is  found  in  the  two-layer  method  (i.e.,  monthly 
rms  difference  and  bias  ^0.7  °C),  the  NPACNFS  shows  less 
accuracy.  As  shown  in  Fig.  11(g)  and  (h),  the  monthly  rms 
difference  is  in  the  range  of  0.8  DC-l,4  °C,  together  with  a 
clear  negative  bias  of  1.0  °C.  From  the  preceding  discussion, 
we  see  that  the  NPACNFS  has  general  underestimation  in  the 
upper  ocean  parameters  (i.e„  D26  and  SST).  This  becomes  a 
critical  factor  to  consider  when  using  NPACNFS  upper  ocean 
parameters  to  study  typhoon-ocean  interaction  since  this  upper 
ocean  is  the  primary  part  of  the  ocean  with  which  typhoons  are 
interacting  [ 25],  [32]. 

Fig.  12  compares  the  profiles  from  the  in-situ,  two- layer, 
and  NPACNFS  profiles  along  a  transect  at  145*  E  for  every  3° 
latitude  from  I0°N  to  25°  N  (i.e.,  passing  through  SNPG  and 
SERZ).  Here,  the  1 45°  E  transect  is  chosen  because  it  is  where 
abundant  in-situ  profiles  exist  with  uniform  distribution  (Fig.  3). 
In  Fig.  12,  one  can  see  that  the  simple  two-layer  profiles  (in  gray 
with  dots)  match  reasonably  well  with  the  observed  profiles 
(black  curves),  while  the  clear  underestimation  ( 1 .0  °C  or  more) 
of  the  upper  ocean  (typically  from  SST  down  to  a  depLh  of 
100  m)  of  the  NPACNFS  profiles  is  seen. 

C  TCHP 

Using  in-situ  profiles.  TCHP  is  calculated  and  used  to  com¬ 
pare  with  the  TCHP  from  the  estimated  profiles  by  the  two- 
layer  and  the  NPACNFS  models.  From  Fig,  13,  one  can  see  that 
the  two-layer  TCHP  has  higher  accuracy  than  the  NPACNFS 
TCHP  throughout  the  three- month  (July-September  2003)  val¬ 
idation  period.  With  the  two-layer  method,  the  rms  difference  is 
about  30  kJ/em2.  and  the  bias  is  less  than  10  kJ/cm2.  In  contrast, 
the  NPACNFS  TCHP  has  evident  underestimation  with  a  large 
negative  bias  ranging  from  -40  to  -60  kJ/cm2  and  the  large 
corresponding  PROM  in  the  range  of  40%-70%  (Fig.  1 3).  This 
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Fig.  1 1  Comparison  0f  the  value  of  the  monthly  mis  difference*  PROM*  and  bias  of  the  NPACNFS  profiles  for  (aHc)  D2Q,  (dHO  P2ft*  and  (g)  and  (h)  SST, 
The  dotted  line  in  lb)  and  (e)  depicts  the  accuracy  threshold  of  the  25%  PROM, 
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V.  Application  in  Studying  the  Intensity  Change: 
of  Supertyphoon  Dianmu  (2004) 

Using  the  TLMJWNP  profiles,  we  study  the  intensity  change 
of  Supertyphoon  Dianmu  (2004).  Dianmu,  which  means 
“Mother  of  lighting”  in  Chinese  and  made  landfall  and  damaged 
Japan,  has  been  the  most  rapidly  intensified  western  North 
Pacific  typhoon  since  I960*2  Within  18  h  (0000  UTC-1800 
UTC  June  15),  it  rapidly  intensified  from  category  I  (70  kn) 
to  the  category  5  superscale  (145  kn).  This  is  equivalent  to 
MK)-kn  intensification  in  24  h*  which  is  more  than  three  limes 
the  usual  “rapid  intensification”  baseline  of  30  kn  in  24  h 
[15],  Observing  the  pre-Dianmu  SSHA  data,  one  can  see  that 
this  happened  a*s  Dianmu  entered  a  large  warm  ocean  feature 
region  in  the  SNPG  that  is  characterized  hy  positive  SSHA 
^20  cm  [Fig,  14(a)].  From  the  TLM_WNP  results,  we  can 
see  that  this  warm  feature  caused  further  deepening  of  D26 
hy  10-30  m  [Fig,  14(b)]  in  the  SNPG,  As  a  result,  it  can  be 


underestimation  is  clearly  the  result  of  the  underestimation 
of  SST  and  D26  in  the  NPACNFS  profiles,  as  discussed  in 
Section  JV-R. 


‘  According  Co  the  l -min  sustained  wind  best  track  data  of  the  U,S.  JTWC 
from  I960  ro  present,  Dianmu  (2004)  and  the  #i  1  typhoon  in  1 98 3  are  1  be  most 
rapidly  intensified  typhoons  in  the  western  North  Pacific;  both  intensified  by 
75  kn  within  18  h. 
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Fig  12.  Comparison  of  ibc  (black  I  in-situ*  (gray  with  dots)  TLM_WNP,  and  (dotted  curve!  NPACNF5  profiles  in  the  SNPG  and  SF:RZ  Along  the  transect  of 
145°  E  from  10-25°  N*  the  profiles  for  every  3°  interval  are  compared.  In  the  13“  N  and  19°  N  TLM_WNP  profiles*  the  SST  data  are  not  available, 


seen  that  this  rapid  intensification  took  place  at  a  region  with  an 
extraordinarily  thick  layer  of  warm  water,  i.e.,  D26  reaching 
150  m  I  Fig,  14(c)]*  The  corresponding  pre typhoon  TCHP  is 
about  120-150  kJ/em2  [Fig*  14(d)],  which  is  higher  than  the 
well-known  Hurricane  Opal  ( 1 998)  case  of  80- 1 20  kJ/em2  [  1 0], 
[32].  As  it  progressed  into  the  interior  of  the  warm  feature 
ai  0600  UTC  June  16,  Dianmu  reached  its  peak  intensity  at 
155  kn  l  Fig,  14(a)  and  (c)].  This  peak  intensity  sustained  for 
12  h  until  1800  UTC  June  16  and  started  to  decline  to  the 
category  4  1 20  kn  at  1 800  UTC  June  1 7  as  Dianmu  moved  into 
a  region  with  a  shallower  warm  layer  of  90  m  [Fig.  14(c)].  This 
situation  can  also  be  seen  from  the  TLM_WNP  profiles*  where 
profiles  from  ibe  intensification  region  (i.e.*  profiles  1-6),  the 
intensity-sustained  region  (profiles  7  and  8),  and  ihc  intensity 
declination  region  (profiles  9-1 2)  are  depicted  in  red,  green*  and 
blue,  respectively  (Fig.  15(a)].  Here,  two  in- situ  GTSFP  profiles 
(profile  6'  and  90  that  were  found  prior  to  Dianmu fs  passing 
(from  June  10  and  12,  2004;  one  near  profile  6  and  one  near 


profile  9)  are  depicted  in  black  (Fig,  15(a),  with  the  location 
depicted  as  stars  in  Fig.  14)*  confirming  acceptable  quality  of 
the  TLM^WNP  profiles*  It  can  also  be  noted  that,  with  the 
two-layer  method,  profiles  along  the  track  can  be  obtained,  but 
with  in -situ  profiles*  such  information  is  not  available.  As  in  the 
Dianmu  case*  only  two  in- situ  profiles  are  found.  As  reported  by 
Lin  et  ai  [21],  this  lacking  of  in-situ  data  along  the  typhoon’s 
tracks  is  common,  particularly,  before  the  deployment  of  Argo 
floats  [11],  [12]  in  2003. 

Thus*  with  the  TLM_WNP  profiles  as  initialization,  together 
with  wind  forcing  from  the  JTWC  best  track  wind*  we  can  run 
an  ocean  mixed-layer  model  [26 j  to  study  the  ocean’s  dynamic 
response  to  typhoon  and  estimate  typhoon-induced  cooling  1 1  ]* 
[3],  [25]*  From  Fig.  15(b)*  one  can  see  the  progressive  increase 
in  mixing  as  the  typhoon  intensifies  (profiles  1-6)  and  reaching 
as  deep  as  120  rn,  as  confirmed  by  the  in-situ  GTSPP  profile 
[black  profile  in  Fig*  15(b)*  with  the  location  depicted  as  a 
triangle  in  Fig,  14]  that  was  acquired  immediately  after  the 
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Eig.  13.  Comparison  of  (he  value  of  the  monthly  rim  difference,  PROM,  and  bias  for  the  TC HP  estimation  using  the  Tt.M_WNP  and  the  NPACNFS  profiles  for 
faHc)  SNPG  and  (dHH  SERZ, 
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Fig.  14  fa)  Composite  of  TOPE  X/Posei  don  and  JasomJ  altimetry  measure- 
menis  for  one  cycle  (ten  days)  between  June  5,  2004  and  June  14,  2004, 
showing  the  pre-Dianmu  SSHA,  (b)  AD26  rTLM_WNP-ikrived  D26  minus 
climatological  D26).  (c)  TE.M_WNP-derivcd  D26,  and  (d)  TCHP  estimated 
from  the  two- layer  method.  The  green  stars  and  black  triangle  depict  the  m-situ 
profiles  that  are  found  before  and  during  Dianmu's  passage.  Typhoon  intensity 
in  the  Saflir-Simpson  scale  is  shown  in  the  legend  of  fa). 
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Fig.  15.  {a)  Initial  profiles  along  Dianmu's  track  that  arc  used  in  the  mixed- 
layer  model  experiment  where  the  tnstiu  profiles,  those  along  Dianmu's 
intensification  location,  those  along  Dianmu's  intensity  maintenance  location, 
and  those  along  Dianmu's  intensity  decay  location  arc  depicted  in  black,  red. 
green,  and  blue,  respectively,  (b)  same  as  (a),  but  from  mixed-layer  model 
except  profile  1 3.  which  obtained  from  ARGO. 


intensification  period  on  June  17,  2004.  From  Fig,  16,  one 
can  also  see  the  corresponding  increase  in  self-induced  cooling 
from  0.6  aC  during  category  [  (at  0000  UTC  June  15)  to 
1.4  °C  at  peak  {i.e„  category  5,  0600  UTC  June  16),  From 
0600  UTC  June  16  onward,  although  Dianmu's  intensity  did 
not  increase,  cooling  still  enhanced  to  2. 1  °C  as  Dianmu  moved 
into  a  region  with  shallower  warm  water  [profiles  9-12  in 
Fig,  15(b)*  Fig.  16],  This  situation  can  also  he  observed  in 
the  SST  maps  from  TRMM  and  AMSR-E  that  were  acquired 
before  [on  June  1 3,  Fig.  1 7(a)]  and  after  [on  June  1 K,  Fig,  1 7(b)) 
Dianmu’s  passing.  From  0600  UTC  June  16  onward,  one  can 
see  clear  enhanced  cooling  to  the  right  to  the  typhoon's  track 
[SST  -  27  *C  in  Fig.  I7(b)|  [I],  [2],  [25],  [36]  as  compared 
with  the  warm  pretyphoon  SST  of  >  29  °C  [Fig.  17(a)].  As 
in  the  earlier  intensification  period  (i.e.,  before  0600  UTC 
June  16),  cooling  is  less  pronounced  [Fig.  17(b)]. 
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Fig.  16.  Relationship  between  five  parameters  (Lc.,  self-induced  SST  cooling,  dynamic  TCHP.  pretyphooti  SST.  pretyphoon  TCHP.  and  vertical  wind  shear)  and 
typhoon  intensity. 


Fig,  1 7.  Composite  of  SST  from  the  TMI  and  AMSR-K  passes  in  (a)  June  1 3. 
2004  for  prc-Diawnu  and  in  (h)  June  18,  2(XW  for  after  Dianmu's  passage. 

Finally,  of  much  interest  for  future  forecast  reference  is 
to  identify  parameters  that  are  sensitive  to  typhoon-intensity 
change.  As  shown  previously,  using  the  derived  profiles  to  the 
mixed- layer  model  self-induced  cooling  during  the  typhoon 
along  the  track  can  be  estimated.  Since  upper  ocean  thermal 
profiles  [Le„  Fig.  15(b)]  can  also  be  obtained,  it  can  he  used 
to  calculate  the  dynamic  TCHP  during  the  typhoon.  Therefore, 
two  parameters  that  include  the  ocean's  dynamic  response,  i,e„ 
the  self- induced  cooling  and  the  dynamic  TCHP  during  the 
typhoon,  can  be  calculated.  It  should  be  noted  that,  without  the 
derived  profiles,  these  two  parameters  could  not  be  calculated, 
as  it  is  impossible  to  have  in-situ  profiles  at  each  point  of 
the  track  (Fig,  14).  Together  with  the  other  parameters,  Le„ 
pretyphoon  SST  (blue  curve  in  Fig,  16  and  the  2-D  map  in 
Fig.  17(a)],  pretyphoon  TCHP  [green  curve  in  Fig.  16,  as  from 
Fig.  14(d) [,  and  atmospheric  vertical  wind  shear1  data  along 
the  track  (brown  curve  in  Fig.  16),  we  compare  the  sensitivity 
of  these  five  parameters  to  typhoon-intensity  change  during 
Dianmu's  drastic  intensity  change  period  (i.e.,  covering  the 

11  Wind  shear  is  esti mated  as  the  difference  between  the  200*  and  850* mb 
European  Centre  for  Medium* Range  Weather  Forecasts  L25(>  resolution  re* 
analysis  wind  data  (hat  arc  averaged  over  a  700  X  700  km  box  that  is  centered 
in  the  typhoon  for  every  6-h  interval. 


rapid-intensification  period  from  category  I  to  5,  followed  by 
the  decay  to  category  4,  from  0000  UTC  June  15  to  1800 
UTC  June  17).  From  Fig,  16,  we  can  see  that  throughout 
the  study  period,  there  is  little  change  in  the  vertical  wind 
shear.  With  the  rapid  intensification  and  decay  of  Dianrnu,  the 
shear  was  typically  about  6  m/s  and  shows  little  variability. 
This  is  also  true  for  the  pretyphoon  SST.  In  comparison,  the 
other  three  parameters,  i.e.,  the  pretyphoon  TCHP,  the  dynamic 
TCHP,  and  the  self- induced  cooling,  have  more  variability. 
As  can  also  be  seen  in  the  relationship  between  the  6-h  in¬ 
tensity  change  and  each  of  the  five  parameters  (Fig.  18),  the 
two  parameters  that  include  the  ocean's  dynamic  responses, 
i.e.,  the  self- induced  cooling  and  the  dynamic  TCHP,  have 
the  highest  correlation  coefficient,  with  R2  0,69  and  0.72, 
respectively  (Fig.  18(a)  and  (b)].  This  suggests  that  there  seems 
to  be  a  relation  between  the  6-h  intensity  change  of  Dianrnu 
during  this  period  and  these  two  parameters.  During  this  study 
period,  the  correlation  coefficient  with  ihe  pretyphoon  SST 
and  the  vertical  wind  shear  is  low,  with  R 2  =0,07  and  0.10, 
respectively  [Fig.  18(c)  and  (e)].  As  for  the  pretyphoon  TCHP. 
R2  —  0.48  shows  more  sensitivity  than  the  pretyphoon  SST  and 
shear  but  is  still  lower  than  the  self-induced  cooling  and  the 
dynamic  TCHP. 

It  should  be  understood  that  typhoon  intensity  is  controlled 
by  many  necessary  atmospheric  and  oceanic  parameters  as  well 
as  the  typhoon's  own  structure  |4],  [35].  It  is  still  a  current 
open  question  in  identifying  sensitive  parameters  in  typhoon 
intensity  and  identifying  the  roles  that  each  parameter  plays 
[4],  [5],  [10],  [35].  The  preceding  sensitivity  comparison  shows 
that  pretyphoon  SST  and  vertical  wind  shear  are  less  sensitive 
than  the  other  parameters.  However,  it  does  not  imply  that 
both  parameters  are  not  important.  The  very  warm  pretyphoon 
SST  in  the  range  of  29  °C“30  °C  and  the  low  wind  shear 
of  about  6  m/s  [8|  are  all  very  favorable  for  typhoon  in¬ 
tensification  but  not  sensitive  to  subsequent  typhoon- intensity 
change  when  the  value  of  the  threshold  of  both  parameters  fi.e., 
pretyphoon  SST  >  26  °C  and  wind  shear  <  10  m/s)  are  satis¬ 
factory.  Unless  the  situation  is  dose  to  threshold  value,  both 
parameters  can  probably  become  sensitive  to  typhoon  intensity. 
In  the  Dianmu's  case,  both  parameters  were  far  beyond  the 
threshold  and  experienced  little  change  (blue  and  brown  curves 
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Fig.  IS.  Intensity  change  of  6  h  of  Dianmu  wjih  the  five  parameters,  t.e.,  fa)  self-induced  SST  cooling,  (b)  dynamic  TCHP.  (c)  pretyphoon  SST,  <d)  pretyphoon 
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in  Fig.  16)  during  the  study  period;  thus,  preiyphw>n  SST  and 
vertical  wind  shear  are  not  as  sensitive  as  the  other  parameters. 
This  study  is  the  beginning  of  investigation  using  two-layer- 
method-derived  profiles  to  systematically  study  these  various 
parameters  in  relation  to  the  intensity  change  of  western  North 
Pacific  typhoons.  With  the  current  analysis,  we  found  that  at 
least  for  the  Dianmu  case,  the  upper  ocean  thermal  structure  is 
highly  critical,  and  the  two  associated  parameters,  i.e.,  the  self- 
induced  cooling  and  the  dynamic  TCHP  during  typhoon,  are 
the  most  sensitive  to  the  6-h  intensity  change  as  compared  with 
all  other  parameters, 

VL  Conclusion 

Using  a  large  set  of  >  5000  in-situ  profiles  from  the 
NOAA/GTSPP  database  during  2002-2005,  this  paper  vali¬ 
dates  the  upper  ocean  thermal  structure  data  that  are  derived 
from  using  satellite  altimetry  and  SST  in  a  two-layer  reduced- 
gravity  ocean  model  that  was  proposed  by  Goni  et  al  |9]  and 
Shay  et  al  [32]  for  the  western  North  Pacific  ocean.  The  vali¬ 
dation  period  is  during  the  May-October  typhoon  season  from 
2002  to  2005,  After  validation,  applicable  zones  and  months 
where  this  method  can  apply  are  identified.  We  found  that  the 
two-layer  method  is  applicable  in  the  southern  and  central  west* 
ern  North  Pacific  (i.e.,  122-170°  E,  18-25°N)  during  most  of 
the  typhoon  months.  However,  it  is  not  applicable  in  the  north¬ 
ern  part  (i.e,,  1 30-170  °E,  25-40°  N).  Further  comparison  with 
the  NRL's  NPACNFS  26-level  hydrodynamic  model  shows  that 
the  NPACFNS  also  experiences  a  similar  regional  dependence. 
In  addition,  dear  underestimation  in  the  upper  ocean  parame¬ 


ters  (i.e.,  SST  and  D26)  is  found.  This  leads  to  significant 
underestimation  in  TCHP  using  the  NPACNFS  profiles  from 
-40  to  -60  kl/cm2,  whereas  using  the  two-layer  method,  the 
TCHP  estimation  reaches  higher  accuracy,  which  is  typically 
about  30  kj/cm2.  Applying  the  two- layer- met  hod -derived  pro¬ 
files  in  studying  the  intensity  change  of  supertyphoon  Dianmu 
(2004)  found  that  the  two  parameters,  i.e.,  the  typhoon's  self- 
induced  cooling  and  the  dynamic  TCHP  during  the  typhoon,  are 
the  most  sensitive  parameters  (R2  ^  0.7)  to  the  6-h  intensity 
change  during  the  study  period.  Therefore,  with  the  availability 
of  the  upper  ocean  thermal  profile  that  is  derived  from  the  two- 
layer  method  in  the  applicable  zones/periods  in  the  western 
North  Pacific,  one  can  conduct  future  analysis  to  systematically 
study  the  control  of  upper  ocean  thermal  structure  in  the  inten¬ 
sity  change  of  western  Pacific  typhoons. 
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